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Abstract 

Switching model with one predator and two prey species is considered. The prey 
species have the ability of group defence. Therefore, the predator will be attracted 
towards that habitat where prey are less in number. The stability analysis is carried 
out for two equilibrium values. The theoretical results are compared with the numer- 
ical results for a set of values. The Hopf bifuracation analysis is done to support the 
stability results. 

1. Introduction 

Among related herbivore species, individuals of smaller species like Dik-dik will be vul- 
nerable to a greater range of predator species and are less likely than larger to be able 
to defend themselves against, or to out run, predators. All small species avoid being 
detected by predators. Smaller species are likely to have to seek carefully their scare, 
scattered food item of high quality and form less cohesive and coordinated feeding groups. 
They live singly or in pairs and find their resources within a defended territory. Because 
they are small and vulnerable, they move and feed cautiously and slowly and never move 
far from cover. The size of the territory is presumably determined by the area that a pair 
can defend and by the availability of suitable food at the season of greatest scarcity. These 
species characteristically remain in one vegetation type in all seasons. 

Larger species individuals feeding upon abundant, evenly dispersed, easily found items, 
are likely to be tolerating low quality food. They form enormous, rather formless, feeding 
aggregations of many thousands of animals. Major predators of zebra, buffalo, kongoni, 
toki and Thomson's gazelle are hyena, wild dog, lion, leopard and cheetas. They form 
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groups for defence against predators and more likely depend upon self-defence, group 
defence, group alertness within a group and speed, to avoid being killed by a predator. 
Dense vegetation and broken terrain disrupt visual communication, and flat open coun- 
try favour it. So such groups are more likely to be found where visual communication is 
favoured and where individuals can conform to the group, speed, and direction of move- 
ment. Unless the group remains cohesive and coordinated, the individual risks becoming 
an outstanding target. At all times individuals in groups must remain in communication 
and their speeds and directions when moving must vary little between individuals. Group 
defence is a term used to describe a phenomenon whereby predation is decreased or even 
prevented altogether by the ability of the prey population to better defend or disguise 
themselves when their number is large. Aggregation tends to reduce the risk of predation 
through simple dilution. Hence, doubling the local density of herbivores while predator 
density remains unchanged could lead to a halving of the mortality risk Hamilton [1], 
Bertram [2]. Of course, this assumes that predators do not seek out areas with very large 
prey density Schaller [3]. Pairs of musk-oxen can be successfully attacked by wolves but 
groups are rarely attacked Tener [4]. There are many examples of group defence - Yang 
and Humphrey [5], May and Robinson [6], Holmes and Bethel [7]. Herds remain well 
coordinated even under attack and individuals may benefit from the alertness and com- 
munication within the herd. Individuals tend to conform with their neighbour activities, 
and many hundreds, even thousands of wildebeest can coordinate rapidly in response to 
an alarm. Large groups also benefit from increased probability of detection of predators. 
The hunting success of lions decline if the group size of prey is large Van Orsdol [8] . Chee- 
tah prefer to hunt single animals. Coursing predators select less healthy, old, sick and 
young prey and those who have lost their herds during migrations due to various reasons. 
Animals in poor condition and without group may reduce vigilance rates. 

Each year, some one million wildebeest migrate across the Serengeti mara ecosystem 
D. Kreulen [9]. The crude cost of this movement, relative to neighbouring resident pop- 
ulations of wildebeest, is a 3% increment in mortality per year Sinclair [10]. The overall 
migratory pattern is thought to be related to food supply, which is itself dependent on an 
uneven distribution of rainfall Grzimek and Grzimek [11], Talbot and Talbot [12], Ander- 
son and Talbot [13], Pennywick [14], Maddock [15]. In the dry season, they need to drink 
Jarman [16] and most home ranges include some riverbank habitat, which is preferentially 
used during that season. At the start of wet season, they move away from the rivers to oc- 
cupy the woodlands where they feed predominantly on new green grass leaves, with some 
herbs, many of these plants may be annuals. As grasses mature and rains cease, their 
diet starts to include progressively more browse and they move into plant communities 
where browse is more. Movements could also be influenced by change in requirements 
for specific nutrients. Kreulen [17] noted that Serengeti wildebeest on their wet season 
range preferred a short-grass over a long grass habitat, and that calcium concentration 
were higher on the short grassland. Migration has also been attributed to the wildebeest's 
dislike of wet and sticky ground Talbot and Talbot [12], Anderson and Talbot [13]. In 
this way herbivores maximize the growth potential of the vegetation through rotational 
grazing where the two concentration areas are sufficiently far apart, the movements are 
called seasonal migrations. 

In predator-prey environment, the predator prefers to feed itself in a habitat for some 
duration and then changes its preference to another habitat. This preferential phenomenon 
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of change of habitat by the predator is called switching. There may be several reasons of 
switching of predators. For example - predator prefers to catch prey species in a habitat 
where they are in abundance but after some duration of heavy predation, when the prey 
species population starts declining, the predator changes its preference to another habitat. 
In this habitat prey species are more in number due to less predation i.e. the predator feeds 
preferentially on the most numerous prey species habitat. This is found to be the case 
when prey species is relatively smaller in size with little or insignificant defence capability 
with respect to predator, like small antelope and cruising predators to locate prey. Many 
examples may be cited where a predator prefers to prey species that is most abundant at 
any time, see Fisher-Piette [18], Lawton et al. [19] and Murdoch [20]. The mathematical 
models which have been generally proposed with such type of predator switching are those 
involving one predator with two prey species, e.g., Holling [21], Takashi [22], May [23], 
Murdoch and Oaten [24], and Raughgarden and Feldman [25], Tanskey [26, 27], Prajneshu 
and Holgate [28], Teramoto et al [29]. 

We are motivated by Serengeti ecosystem which can be divided into two habitats for 
wild life - (a) open southern grasslands with low rainfall that support a relatively low 
biomass of short-growing grasses and (b) wooded northern grasslands with higher rainfall 
that support tall, highly lignified grasses Braun [30], McNaughton [31, 32]. Rainfall is a 
key factor influencing primary productivity in both grassland types Braun [30] Sinclair 
[33]. McNaughton [31, 32] . All of the migratory species (wildebeest, zebra, and Thom- 
son's gazelle) show similar seasonal shifts in habitats, using short grassland in the south 
during wet season and tall grasslands in the north during dry season Pennywick [14] and 
Maddock [15]. Their long range movements roughly correspond to seasonal transitions. 
We have considered a system having a predator species interacting with same prey species 
in two habitats. Prey is able to migrate among two different habitats at some cost to the 
population in the sense that the probability of survival during a change of habitat may 
be less than one. The predator can feed on either habitats. The prey species in both 
habitats have the ability of group defence but it will be effective in the habitat where the 
population of prey is large. The predator will be attracted toward s that habitat where 
prey are less in number. 

Freedman and Walkowicz [34] considered a predator-prey system in which the prey 
population exhibits group defence. They described that if the carrying capacity of the 
prey population is sufficiently large and there is no mutual interference among predators 
then the predator population always goes to extinction. Freedman and Quan [35] studied 
predator-prey system with third population extinction. Shigui Ruan and Freedman [36] 
analyzed group defence in Gause-type models for three species. They gave criteria for 
persistence when there is no mutual interference and when there is mutual interference 
of predators. Freedman and Shigui Ruan [37] have discussed a three species food chain 
model with group defence. They have shown that model undergoes a sequence of Hopf 
bifurcations, using the carrying capacity of the environment as a bifurcation parameter. 

This paper is organized as follows - Model formulation is in Section 2 and equilibrium 
and stability is discussed in Section 3. Section 4 includes Hopf bifurcation analysis with 
respect to 5 (conversion rates of prey to predator) and u (death rate of the predator). 
Bifurcation points and effect on stability for selected numerical data values are in Section 5. 
Final discussion and results are summarized in Section 6. 
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2. Prey-Predator model with group defence 

We consider a class of Voltera predator prey model with group defence exhibition by the 
the prey of the form 

*i(0)>0, 

X 2 (0) > 0, (2.1) 
^(0) > 0, 

where ^- = —J , and X±(t), X 2 (t) are prey species in the first and second habitats re- 
spectively and Y(t) denote predator species which feed upon X\ and X 2 . Prey species are 
able to disperse among 2 different habitats at some cost to the population. Prey have the 
ability of group defence, so predator will go towards the habitats where prey population 
is less numerically. We consider that in the beginning, prey species are less in habitat 
1 and so predator will be attracted towards them. For giving protection to their mates 
the prey species of the habitat 2 will rush towards habitat 1. In this way their strength 
in the habitat 2 will fall short and the predators will be attracted towards them. So to 
protect them prey start coming from hibitats 1 to 2. This implies a kind of switching from 
one habitat source of food to another habitat as the prey fall short alternately. Due to 
seasonal migration of prey species, none of the habitat population will be extinct. This 
situation is described by the above model where 

fli : Predator response rates towards the two prey X\ and X 2 respectively. 

5i : The rate of conversion of prey to predator. 

€i : Inverse barrier strength in going out of the first habitat and the second habitat. 
Pij : The probability of successful transition from i th habitat to j th habitat. 
cti : Specific growth rate of the prey in the absence of predation. 
v : Per capita death rate of predator. 

We assume that Pi,5i, 6i,pij, cti, and v are positive constants and e\p\ 2 > /3 2 Y(0), 
C2P21 > PiY(0), so that the X\ and X 2 are not negative. 

3. Analytical Solution 

In this section, we proceed to analyze the system (2.1). We examine the equilibrium of 
this system. We obtain equilibrium solutions by setting time derivatives to zero. If in 
equations (2.1) Xi,X 2 ,Y are the equilibrium values of Xi,X 2 ,Y respectively, there are 
two possible equilibria, namely, 

(i) X\ = X 2 = Y = 0, i.e., population is extinct and this always exists. 



[3\X Y 

Xi = (ai - ei)Xi + e 2 p2iX 2 - 2 ' 



X-2 = (a 2 — e 2 )X 2 + e\p\ 2 X\ 



(Xi + x 2 ) 

(3 2 XfY 
(X! + X 2 ) 



• 1 u+ (X 1 + X 2 ) + (X 1 + X 2 ) ] ' 
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(ii) x 1=f + x 2 = "V + V 



(SiPi + S2P2X 2 ) ' (<5i/3i + 5 2 p 2 X 2 ) ' 

v (1 + X) {(a 2 - e 2 ) + e lPl2 X} . 

y = — — or equivalently 

(1 + X) {( ai - ei )X + e 2P21 } 

Pi 

Here, X = (X1/X2) is a real positive root of the cubic equation, 

2 ( ai - ei)X 3 + /3 2 e 2 p 21 X 2 - faewsX - /?i(a 2 - e 2 ) = 0. (3.1) 

For equilibrium values (Xi,X 2 ,Y) to be positive, a positive real root of (3.1) must be 
bounded as 

< X < (3.2) 
eiPi2 ei - «i 

Since in the absence of the predator — > and — — ^ > 0. 

at at 

From Descarte's sign rule the equation (3.1) has unique positive root and hence the 
unique positive equilibrium solution of model (2.1). 

(a) Stability Analysis of Equilibirum (i). Consider a small perturbation about the 
equilibrium level X\ = X\ + u, X 2 = X 2 + v, Y = Y + w. Substituting these into the 
differentiate equation (2.1) and neglecting products of small quantities, we obtain stability 
matrix 



(a\- e\ e 2 p 2 i \ 

£1^12 a 2 — e 2 
\ -v J 



(3.3) 



The characteristic equation of this matrix is 

[y + a) [A 2 - A{qi - ei) + (a 2 - e 2 )} + {cx\ - ei)(a 2 - e 2 ) - eie 2 pi 2 p 2 i] = 0. (3.4) 

If a\ < e\ and a 2 < e 2 and (ai — e\)(a 2 — e 2 ) > e\e 2 pi 2 p 2 ±, the equilibrium X\ = X 2 = 
Y = is locally stable otherwise unstable. 

(b) Stability Analysis of Equilibrium (ii). Following the above procedure as in (a) 
the stability matrix becomes: 

( A -Si A \ 



X 2 {X l + X 2 ) 

x 2 _ B B p 2 xl 



X! ' (X! + X 2 ) 

\ CY DY 



(3.5) 
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The characteristic equation associated with the positive equilibrium (X±,X 2 ,Y) of this 
model is 



A 



X{A + B) + X {WTx 2 ) + WTx 2 j) 



( AC(3 2 YX\ BCfcYXl ADY (3 2 X\ BDY^Xj ., 
\X 2 (X 1 + X 2 ) (X l +X 2 ) {X 1 +X 2 ) X 1 {X 1 +X 2 ) ] 



where 



(hXlY (3 2 XfY 
A = (cci - ei) + . - ■ % . , J B = (a 2 -e 2 ) + 



r .,m CftYX^ DYfoX* 
ai = -(A + B), a 2 = —= + 



(3.6) 



(X1+X2) 2 ' V y (^i+^ 2 ) 2 ' 

_ (_ y + 2Xij 2 /3 2 )y z? _ (-1/ + 2x2^1 )y 

(*l+* 2 ) ' X!+X 2 ) 

Equation (3.6) can be written in the form 

A 3 + aiA 2 + a 2 A + a 3 = 0, (3.7) 

where 



(3.8) 



(X, + X 2 ) (X, + X 2 ) ' 
/ ACfcYXf BCpiYXl ADY(3 2 Xf BDYfi\X\ 

° 3 " " \x 2 {x x +x 2 ) + (xl + x 2 ) + (Xi + x7) + Xi(Xi + x 2 ) 

The Routh-Hurwitz stability criteria for the third order system is 

(a) 01 > 0, a 3 > 0, 

(b) a\a 2 > a 3 . 

Hence, the equilibrium (ii) will be locally stable to small perturbations if it satisfies the 
following conditions 

A + B < 0, 

(CXt + DX 2 ) {A(3 2 Xt + BfaXl) < 0, (3.9) 
(/?iXf - foXf) (BDX 2 - ACX X ) > 0. 

Here we observe that the stability of the model depends upon the conditions (3.2) and 
(3.9) together with various parameters. 



4. Hopf Bifurcation Analysis 

We investigate the Hopf bifurcation for the system (2.1) taking 6\,5 2 and v as the bifur- 
cation parameters. 

First, we determine the criteria for Hopf bifurcation using <5i (rate of conversion of prey 
in the habitat 1 to predator) as the bifurcation parameter. For non zero equilibrium we 
look at the characteristic equation (3.7). 



For A < 0, B < 0, C > 0, and D > 0; 
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ai, d2 and as are positive in (3.8), clearly (3.7) has two purely imaginary roots iff aia 2 = as 
for some value of 5± (say 5± = Si). There exists a unique Si such that aia 2 = as- 
Therefore there is only one value of Si at which we have a Hopf bifurcation. Thus in the 
neighbourhood of Si the characteristic equation (3.7) cannot have real positive roots. For 
Si = Si, we have 

(A 2 + a 2 ) (A + ai) = 0, (4.1) 
which has three roots 

Ai = iy/ai, A 2 = -i\fa^, A 3 = -ai. 
The roots are in general of the form 

\i(Si)=p(Si)+iq(Si), X 2 (Si)=p(Si)-iq(Si), Xs(Si) = -ai(Si). (4.2) 

To apply Hopf's bifurcation Theorem [38] to (2.1), we have to verify the transverality 
condition 

Substituting Xk(Si) = p(Si) + iq(Si) into (3.7) and calculating the derivative, we get 

R(S 1 ) P '(S 1 )-S(S 1 )q'(Si)+T(Si) = 0, 
S(Si) P '(Si) + R(5i)q'(Si) + U(5i) = 0, 

where 



(4.4) 



R(Si) = Sp 2 (Si) + 2ai{Si)p{Si) + 02(^1) - 3q 2 (S 
S(Si) =6p(Si)q(5i) + 2ai(Si)q(Si), 

T(Si) = p 2 (Si)a[(Si) + a' 2 (Si)p(Si) + a' 3 (Si) - a'^q^Si), 
U(Si) = 2p{Si)q{Si)a' 1 (Si) + a' 2 (Si)q(Si). 



(4.5) 



If SU + RT / at Si = Si, then 

'd\ k \ (SU + RT) 



Re 



+ 0. (4.6) 

<5i=5i 



dSiJ Sl= - Sl 2(Ri + Si) 
Now from equation (4.5) 

SU + RT = ai a 2 - a's at S 1 = S 1 , (4.7) 

where a' 2 = -77^ and a's = -rr~- X is a real positive root of equation (3.1) and independent 
ddi ddi 

oiSi ^ = 0. 
dSi 

Substituting the values of 01, a' 2 and 03 in equation (4.7), and using aia 2 = as at Si = Si 
which gives 

( pa2 ~ §0 (it " AC ) = at 61 = Sl > 
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i.e., -— = X, at Si = 5± and value of 1 , we obtain, 
AC do\ 

<sc, + flT) = (Tw(|-*^) BAi ' 2 

V _ _ 7 _ (4.8 

*2<?i/?i -v + X 1 5 2 P 2 + X(-v + 2X X S 2 (3 2 

-v + 2X x b 2 fo 

Since C and D are positive, the terms in square bracket are positive. Hence 

Re(§) if ^ + X\ k=l,2. 

\d5iJ Sl= - Sl (3 2 

and A 3 ((5i) = -ai(<Ji) + 0. 

We summarize the above results in the following: 

Theorem. Suppose E(X 1 ,X 2 ,Y) exists and A < 0, B < 0, C > and D > 0. TTie 
system (2.1) exhibits a Hopf's bifurcation in the first octant leading to a family of periodic 

Pi -s 

solutions that bifurcates E for a suitable value of 5\ in a neighbourhood of Si if — / X . 

P2 

We can get a similar result when S 2 is taken as a bifurcation parameter. Therefore the 
bifurcation points that we obtain in Table I of section 5 are the Hopf's bifurcation points. 

Now we analyze the dynamics of (2.1) with respect to v (per capita death rate of 
predator). 

Similar to (4.3) we need to verify 



Re(^] SU + RT /0 fc-12 



(4.9) 



where S, U, R, T have the similar expression as given in (4.4) but now these are functions 
of v (instead of <5i). 

SU + RT = aia 2 — a' 3 , where a 2 = and 03 = 

Using the relation 

dXi (X + 1)X 



(4.10) 



du (Si/3i + S 2 p 2 X 2 Y 
we can show that 

For purely imaginary roots of (3.7), aia 2 = 03 at v = v, we get 

{A + B) + D(3 2 X^J = AC[3 2 X 2 + + ADX. (4.12) 

Now using (4.10)-(4.12) we obtain the value of 

aia' 2 - a' 3 = 0. (4.13) 
Therefore, from (4.9), we get 

Re ijt) -=°' fc = 1 ' 2 ' 
i.e., if v is a bifurcation parameter, there is no Hopf's bifurcation. Hence with respect v 
the system (2.1) is either stable or unstable, the result which we obtain in Section 2. 
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5. Numerical Solutions 

For illustration we have seen the effects of /3's, c)'s. The behaviour of stability with respect 
to /3's and <5's is given in Table I. 

Table I 



0.01, 


CKl = 0.05, Ct 2 = 


0.25, ei = 


= 0.1, e 2 = 


0.3, P12 = 0.5, P21 


01 


02 


Si/S 2 


Bifurcation point 


stable 


0.01 


0.02 


5 2 = 0.3 


Si = 


0.75123 


0.75124 < 61 < 1 


0.01 


0.02 


5 L = 0.5 


62 = 


0.19968 


< 8 2 < 0.19967 


0.01 


0.03 


<5i = 0.5 


62 = 


0.12949 


< 8 2 < 0.12948 


0.02 


0.01 


5 2 = 0.3 


Si = 


0.16448 


< 5i < 0.16447 


0.01 


0.01 


Si = 0.5 


82 = 


0.42507 


< 8 2 < 0.42506 


0.015 


0.01 


5 L = 0.5 


62 = 


0.66456 


0.66457 < 8 2 < 1 


0.02 


0.01 


Si = 0.5 


82 = 


0.91196 


0.91197 < 5 2 < 1 



The last 2 columns give the values of 61/82 at which the model is stable or unstable. 
In fact fourth column gives the bifurcation point (where the model is stable, below/above 
that value the model is unstable/stable). In the appendix we show that these bifurcation 
points are in fact Hopf bifurcation points. We have also done computations to see the 
effect of v on the stability when 0\ > 2 . It is interesting to find that the model is either 
stable or unstable with respect to 5's. and Hopf bifurcation does not exist with respect 
to v. 

One can try similar analysis with respect to other parameters also. 
The set of equations given in (2.1) have been numerically integrated for four cases given 
in Table II with other parameters as given in Table I. 

Table II 



Case 


0i 


02 


Si 


82 


Stable/unstable 


(i) 


0.01 


0.02 


0.95 


0.3 


stable 


(ii) 


0.02 


0.01 


0.5 


0.3 


stable 


(iii) 


0.01 


0.02 


0.15 


0.3 


unstable 


(iv) 


0.02 


0.01 


0.6 


0.3 


unstable 



These sets were picked up while doing the computations of analytical results in previous 
section where the behaviour of the model is shown in the last column. The initial condi- 
tions used are the corresponding equilibrium values in each case with slight perturbations. 
Figures 1 to 4 give the behaviour of Xi,X 2 and Y with respect to t in above four cases 
and as expected we get stable behaviour in the figures 1 and 2 and unstable behaviour 
in the figures 3 and 4. Figures 5 to 8 give the prey-predator dynamics when the model 
is stable whereas figures 9 to 12 represent the prey-predator dynamics when the model is 
unstable. These contours also support the predictions of Table II. 
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6. Summary and Discussion 

We have considered a system having a predator species interacting with prey species in 
two habitats. Prey is of large size and migrate between two different habitats at some cost 
to its population in the sense that the probability of survival during a change of habitat 
is less than one. The predator can feed on either habitats. The prey species in both 
habitats have the ability of group defence but it will be effective in the habitat where the 
population of prey is large. Due to group defence ability of the prey, predator will select 
the habitat where prey might have insufficient defending capability (i.e. numerically less, 
old, sick, some youngs and those who might have lost their group during migration due 
to various reasons). The stability analysis has been carried out for both zero and nonzero 
equilibrium values. Nonzero equilibrium X\ and X2 for prey in first and second habitats 
depends on the death rate of predators i.e. if the death rate of the predator is high then 
values of X\ and X2 will increase or vice versa. X\ and X2 values will decrease if the 
predator response towards both habitats increases respectively. Nonzero equilibrium will 
be stable if it satisfies all three conditions of equation (3.9). Consider first the limit where 
ei — > and 62 — ► 0, there is no movement of prey in or out of both habitat and so values of 
A and B in (3.9) 1; will be positive and nonzero equilibrium will become unstable. Hence, 
increasing the values of ei and €2 always stabilizes. We can conclude from (3.9) 2 that 
stability increases for increasing values of conversion rate of prey species by predator. 
From (3.9)3 we see that if X2 > X\ then for stability it is more likely (3\ > P2 provided 
BD 

> X or vice versa i.e. if the equilibrium value of prey species in second habitat 

is more than first habitat then predator will attract towards first habitat because prey 
exhibits group defence and our model predicts this behaviour. 

Hopf bifurcation analysis has been carried out for both models with respect to A (con- 
version rate of prey to predator) as a parameter. In the sense of ecology, Hopf bifurcation 
has helped us in finding the existence of a region of instability in the neighbourhood of 
nonzero equilibrium, where prey species in both the habitats and predator will survive 
undergoing regular fluctuations. However, the conditions of Hopf bifurcation might not 
be satisfied due to changes in other parameters and change the steady state or otherwise. 
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